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GRAPHICAL METHOD FOR PREDICTING LIFE OF A ROCKET THRUST CHAMBER 
WITH HALF-HARD ZIRCONIUM COPPER LINER AND 
ELECTROFORMED NICKEL CLOSEOUT 
by Harold J. Kasper 
Lewis Research Center 

SUMMARY 

A simple method for estimating the life of a regeneratively cooled rocket thrust 
chamber has been developed and is based on the hot -gas wall temperature and the tem- 
perature difference between the hot-gas wall and the outside surface of the closeout. 

This method permits a quick estimate of the life of a thrust chamber when design 
changes or test-cycle variations are considered. Strain range and life are presented 
graphically as functions of these temperature parameters for a typical high-performance 
rocket thrust chamber with a half -hard zirconium copper liner and an electroformed 
nickel closeout. 


INTRODUCTION 

The Space Shuttle has introduced many new design requirements for rocket thrust 
chambers. One of these is reusability, which has not been a major design requirement 
in the past. Now that more emphasis is being put on economy in space endeavors, the 
need for reusability is apparent. In order to minimize engine weight and volume re- 
quirements, future rocket engines may be required to operate at chamber pressures well 
above those of present engines. Consequently, the heat flux in the throat region will in- 
crease considerably. 

One of the chief concerns in the design of a reusable regeneratively cooled thrust 
chamber is the fatigue life of the thrust chamber liner. The increased heat flux in the 
throat region, which is a result of higher chamber pressures, contributes significantly 
to the complexity of this design problem. The large temperature differences that occur 
between the inner wall and the outer shell during engine operation usually create plastic 
strains that have a major effect on thrust chamber life. During engine startup when the 



inner wall is hotter than the outer shell, large compressive plastic strains are created 
near the inner wall. When the engine is shut down and if the coolant flow continues in the 
cooling channels, which is typical for a fatigue-life test cycle, the inner wall becomes 
cooler than the outside shell and tensile strains occur near the inner wall. The range 
between these strains is referred to as the total strain range. There is an inverse rela- 
tionship between this total strain range and cyclic life: as total strain range increases, 
cyclic life decreases. Therefore, in order to design an engine for a predetermined life, 
it is necessary to predict the cyclic strains as accurately as possible. 

A program has been initiated at the Lewis Research Center to study the effects of 
thrust chamber liner materials, coatings, and designs on chamber life. Low strength, 
high conductivity materials in particular are being studied because the theoretical 
coolant wall thickness required by the higher strength, low conductivity materials to 
transfer the increased heat load at acceptable material temperatures is too small to sus- 
tain the coolant pressure load. Thrust chambers employing the various materials are 
fabricated and then subjected to a cyclic combustion environment until failure occurs. In 
conjunction with this, the cyclic lives of these thrust chambers are predicted so that the 
analytical methods can be verified. During the course of the test program, it is desir- 
able to have a simple method of predicting life based on measured wall temperatures so 
that if the test cycle conditions are adjusted, the corresponding effect on life can be 
quickly estimated. 

The purpose of this report is to present a simple graphical method of predicting the 
life of a typical test thrust chamber based on the hot -gas wall temperature and the tem- 
perature difference between the hot -gas wall and the closeout outside surface. The 
graphs are for a test thrust chamber with a 6. 6-centimeter (2. 6-in. ) diameter throat, a 
half -hard zirconium copper liner, and an electroformed nickel closeout. The propellants 
used are hydrogen and oxygen at an oxidant-fuel ratio of 6. 0 and a chamber pressure of 
4. 137 meganewtons per square meter (600 psia). A more detailed description of the 
chamber is given in reference 1. The method of analysis and the resulting graphs are 
presented herein. 


ANALYTICAL METHOD 

The conventional analytical approach to life prediction of a regeneratively cooled 
thrust chamber such as the one shown in figure 1 includes a computer thermal analysis 
and a computer structural analysis of the throat section. The throat section is normally 
chosen for the life analysis because the high heat flux that occurs here generally pro- 
duces the maximum strain. Temperatures obtained from the thermal analysis program 
are input to the structural analysis program as loading conditions along with the thrust 
load and the pressure loads. 
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The thermal map required for the structural analysis is usually generated with a 
thermal analyzer computer program such as SDSfDA (ref. 2). Input to this program re- 
quires hot -gas and coolant wall heat -transfer boundary conditions as functions of time. 
Temperatures such as those in figure 2 are output for the selected cross section at arbi- 
trarily selected time intervals. The numbers appearing on the model in figure 2 identify 
the locations at which temperatures are calculated, and the accompanying table lists 
typical corresponding temperatures. Figure 3 shows the throat section structural model 
and the identifying element numbers for the structural finite -element analysis. The tem- 
perature map required by the structural analysis program is based on the structural 
model; that is, a temperature is required at the center of each structural element of fig- 
ure 3. Figure 2 shows that the SINDA thermal map provides 68 internal temperatures 
and 23 boundary temperatures. Therefore, the SINDA temperature that is nearest to the 
center of each of the 34 structural elements is selected as the temperature of that el- 
ement. In some instances, two temperatures near an element center are averaged. A 
revised temperature map corresponding to the structural element center locations is 
shown in figure 4. 

Although the SINDA program provides temperatures as a function of time, only two 
times are selected for the structural analysis. These times correspond to the points of 
maximum temperature difference AT between the hot -gas wall and the closeout outside 
surface during the start transient and during the shutdown transient of the firing cycle. 
Generally, this is when the maximum compressive and tensile strains occur. Diming the 
start transient, large compressive plastic strains are induced in the copper liner. This 
is a result of the hot copper liner being restricted in its thermal growth by the cooler 
nickel closeout. The magnitude of these strains is related to AT. During the test cycle 
shutdown transient, tensile strains are induced in the copper liner because the coolant is 
allowed to flow continuously as the engine is started and stopped, which allows the close- 
out to acquire a higher temperature than the liner. The relationship between coolant 
flow, cyclic AT and chamber pressure is shown in the two cycle schematics of figure 5. 
Both firing cycles were used during the test program and each produced approximately 
the same temperature differences during the start and shutdown transients. 

Once the thermal maps are obtained, the fatigue-life strain range of each element is 
calculated by using the thermal map for the start transient peak AT as input to the 
three-dimensional finite -element program RETSCP (refs. 4 and 5). Total strain and re- 
sidual strain (plastic strain) are determined for each element. Residual strains are then 
input along with the shutdown transient peak AT thermal map, and total strains for the 
combined loading are provided as output. The strain range for these loading conditions 
is obtained from a typical stress -strain hysteresis loop (fig. 6). During the start tran- 
sient, the total strain in the critical element of the hot liner includes plastic de- 
formation. During the shutdown transient, the state of strain in the critical element re- 
verses along the elastic -plastic line until a total strain of e£ ot is achieved. In the 
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process of this reversal the strain state passes through a zero stress point that corre- 

XJ 

sponds with the residual plastic strain e p . The strain range is the strain between the 
peak strain points of the loop and is given by 

Strain range = Ae = ef Qt + eg* - e p (1) 

Thrust chamber life is then predicted from a curve that relates strain range and 
cycles to fracture for the liner material. Such a curve for half -hard zirconium copper 
was taken from reference 6 and is shown in figure 7. 


GRAPHICAL METHOD 

The graphical method eliminates the thermal and structural computer analyses and 
allows the strain range and hence life to be determined from the hot -gas wall temper- 
ature and the temperature difference AT. The method used to construct these graphs 
and their application follows. 


Thermal Map Development 


In order to produce a family of curves that show the relationship between strain, hot- 
gas wall temperature, and temperature difference, it is necessary to compute a thermal 
map for each point needed to plot a curve. This becomes a lengthy procedure when using 
SINDA because input boundary conditions must be repeatedly adjusted until the desired 
temperatures are output at the hot-gas wall and the closeout outside surface. Therefore, 
temperature equations and weighting factors were used in place of SINDA to generate the 
required thermal maps. 

The temperature equation for the critical time ’during the start transient was taken 
from reference 3 and modified. This resulted in 


where 



T temperature at radius r 


( 2 ) 
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Tj hot -gas wall temperature 
T q closeout outside surface temperature 
r Q outside radius 
r variable radius 

inside radius 

This equation provides a radial temperature distribution only and does not account for 
any circumferential variation. Therefore, the temperature calculated with this equation 
for each row of structural elements was considered as an average temperature for that 
row. 

To account for the circumferential variation in temperature, a weighting factor was 
applied to the temperature calculated from equation (2). The weighting factors were de- 
termined from SINDA thermal maps by calculating the ratio of each element temperature 
to the corresponding average temperature for that element row. For example, referring 
to figure 4, the average temperature for the element 30 to element 34 row is 728 K 
(1311° R). The weighting factor for element 30 is the ratio of 749 K (1348° R) to 728 K 
(1311° R) or 1. 028. This was done for each structural element of the revised SINDA 
thermal maps associated with the cases of table I. (These cases were analyzed during 
the development of the final test cycle. ) Hence, there were five weighting factors for 
each element based on the five thermal maps. The final weighting factor for each el- 
ement was taken as the average of these five factors. 

In examining the SINDA thermal maps, it was observed that the circumferential tem- 
perature variation was approximately the same for the corresponding rows of the differ- 
ent cases. Therefore, it was felt that the average of the circumferential weighting fac- 
tors would adequately account for the circumferential temperature variation of all cases. 
Furthermore, the radial temperature gradient is usually predominant in influencing the 
thermal strain throughout the cross section. The weighting factors and the average 
weighting factor for each element are given in table n. 

To calculate a thermal map for a particular hot -gas wall temperature and temper- 
ature difference, the midpoint temperature for the element row was calculated with equa- 
tion (2). This value was multiplied by the weighting factor for each element in that row 
to obtain the element temperature. This procedure was executed for the various element 
rows appearing in the structural model. 

The procedure used to generate the maximum strain thermal map for the shutdown 
transient was the same except for the radial temperature equation. Two equations were 
used in place of equation (2). They were 

T = Tj (3) 
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where 3. 30 s r < 3. 38 centimeters (1. 300 ^ r < 1. 3325 in. ), and 


T = (T. - T q ) 



+ T, 


( 4 ) 


where 3. 38 < r <4.18 centimeters (1. 3325 < r < 1. 6440 in. ) and r- = 3. 38 centimeters 
(1.3325 in.). 

The temperature variation between the hot-gas surface (r = 3. 30 cm; 1. 300 in. ) and 
the coolant surface (r = 3. 38 cm; 1. 3325 in. ) at the time of the shutdown maximum tem- 
perature difference was small. Therefore, all element temperatures within these bound- 
aries were considered to be equal to the hot-gas wall temperature (eq. (3)), and the in- 
side radius of equation (4) was taken as 3. 38 centimeter (1. 3325 in. ). The weighting fac- 
tors for the shutdown transient temperatures are tabulated in table ID. A discussion of 
the accuracy of the equation-factor method for calculating thermal maps is included in 
the appendix. 


Graph Construction 

Structural computer runs were made using start transient thermal maps derived 
from the equation -factor method with various inside and outside temperatures. Max- 
imum strain was plotted as a function of AT and hot -gas wall temperature in figure 8. 
Notice that the upper limit of each curve occurs at a AT equal to the hot-gas wall tem- 
perature, which corresponds to an outside temperature of 0 K (0° R). Separate shutdown 
transient structural computer runs were made independent of the start transient runs us- 
ing equation calculated thermal maps and neglecting residual strains. Maximum strains 
for the shutdown transient runs are plotted in figure 9. The AT values along the ab- 
scissa of the graph in this figure are shown as negative values. This is because the 
chamber outside surface temperature is greater than the hot-gas wall temperature during 
shutdown and because the previous definition of AT has been maintained (AT = T^ - T Q ). 
To estimate the strain range, the strain corresponding to the start transient AT and the 
hot -gas wall temperature is read from figure 8 and the strain corresponding to the shut- 
down transient AT and the hot -gas wall temperature is read from figure 9. The two 
strains are added to obtain the estimated strain range. 

The graph in figure 10, which relates cycle life to AT and hot-gas wall temper- 
ature, is a composite of figures 7 to 9. Strain values corresponding to various AT 
values and hot-gas wall temperatures were taken from figure 8 and increased by 
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10 percent to account for the small shutdown transient strain that is added to the start 
transient strain to obtain the strain range. The shutdown transient strain is normally 
considerably smaller than the start transient strain because the maximum shutdown 
transient AT is much less severe. Cycle lives corresponding to the increased start 
transient strains were taken from figure 7 and plotted against the corresponding AT 
values and along corresponding constant hot -gas wall temperature lines. Hence, an esti- 
mated cycle life can be determined directly from AT and the hot -gas wall temperature. 


Application of Graphs 

As an example of the application of the graphs, a start transient hot -gas wall tem- 
perature of 833 K (1500° R) and a AT of 667 K (1200° R) are assumed. A maximum 
start transient strain of 2. 38 percent is obtained from figure 8. A maximum shutdown 
transient strain of 0.25 percent, corresponding to an assumed hot -gas wall temperature 
of 111 K (200° R) and a AT of -167° C (-300° F) is obtained from figure 9. The strain 
range is then calculated by adding these two strains, which results in a strain range of 
2. 63 percent. Using this value and figure 7 results in an expected fatigue life of 770 
cycles. An expected fatigue life of 770 cycles is also obtained from figure 10 with the 
same start transient temperatures that were assumed previously. 


RESULTS AND DISCUSSION 

I 

In an effort to evaluate the effect of neglecting the residual strains in the graphical 
method, a structural analysis that included the residual strains was performed. The 
thermal map corresponding to the start transient AT and hot-gas wall temperature of 
the example given in the previous section was input to the structural program. The re- 
sidual strains from this analysis were then input along with the shutdown transient AT 
and hot -gas wall temperature of the same example. The resulting strain range was 2. 44 
percent, while that obtained from the graphical method was 2. 63 percent. Hence, the 
inclusion of residual strains had a little effect on the total strain range. Furthermore, 
the expected fatigue life for the 2.44 percent strain range was 900 cycles and for the 
2. 63 percent strain range, 770 cycles. Therefore, the graphical method is conservative 
because the larger strain range obtained by adding the graphical values resulted in a 
shorter expected life. 

Table I lists the cycle conditions of five theoretical firing cycles used in the develop- 
ment of the test program firing cycle. Table IV shows maximum strain range values for 
these five cases. The values in the second column were obtained with SINDA generated 
thermal maps and RETSCP and have residual strains included. Also shown are 
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maximum strain range values obtained from figures 8 and 9. Although the particular lo- 
cation of the maximum strain range varied slightly, it was always within the 0. 09 centi- 
meter (0. 035 in. ) inside annulus of the chamber where the strain ranges were generally 
much larger than those in other locations. This typical area of large strain ranges can 
be seen in the strain range distribution of case 5 (fig. 11). The strain ranges obtained 
from the graphs are larger than those with the residual strains included in four of the 
five cases as seen in table IV. It was again concluded from this and the fact that life 
decreases with increasing strain range, that the graphical method will normally give 
more conservative life predictions than those obtained with the more rigorous method. 

Figure 10 shows that cycle life decreases with increasing AT and increasing gas 
wall temperature. Also, at the lower AT values (below approximately 555 K (1000° R)) 
and the lower hot-gas wall temperatures (below approximately 722 K; 1300° R), it ap- 
pears that an equal reduction in either of these has about the same effect on life. At the 
higher AT values and the higher hot -gas wall temperatures, the slopes of the curves 
decrease and a reduction in AT seems to have less of an effect on increasing life than 
an equivalent decrease in the hot -gas wall temperature. This is probably due to the ma- 
terial life curve. At the lower strain ranges (below approximately 2%), which corre- 
sponds to the lower AT values and the lower hot -gas wall temperatures, cycle life in- 
creases more rapidly than above 2 percent for an equivalent change in strain range. 


CONCLUDING REMARKS 

A simple method for predicting the life of a regeneratively cooled thrust chamber 
with a half -hard zirconium copper liner and a nickel closeout has been presented. 

Graphs that relate strain range and life to the hot-gas wall temperature and the temper- 
ature difference between the hot-gas wall and the closeout outside surface were prepared 
by substituting temperature equations for the computer analyzer and by neglecting the re- 
sidual strains in the structural analysis. The accuracy of the graphs was checked by 
comparing results obtained from the graphs with those obtained from the more sophis- 
ticated computer analysis, which employed SINDA and included residual strains in the 
structural analysis. Since the results from the two methods were fairly close, it may 
be concluded that figures 7 to 10 are satisfactory for quick strain range and life esti- 
mates of thrust chambers similar to those being tested at Lewis. 

Although the simplified analysis is for a half -hard zirconium copper liner and a 
nickel closeout, it should be possible to estimate strain ranges and lives of other liner 
materials with nickel closeouts if the thermal coefficient of expansion is approximately 
equal to that of half -hard zirconium copper. Mechanical properties, such as yield 
strength and the modulus of elasticity effect the total strain range; however, the thermal 
coefficient of expansion is usually the most influential. The mechanical properties have 
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more of an effect on the fatigue life curve so that the appropriate life curve must be used 
once the strain range is estimated from figures 8 and 9. 

Variations in chamber pressure and coolant inlet pressure also effect the strain 
range and, hence, life. Small variations of 10 to 15 percent in either pressure should 
have only a small effect on the total strain range when compared with the effect of the 
thermally induced mechanical strain at high temperatures. As the hot-gas wall tem- 
perature approaches the melting point of the liner material, coolant pressure effects 
become more predominant because the liner material strength properties are reduced 
with increasing temperature. Furthermore, as the temperature increases, failure 
modes other than fatigue could become the primary failure mechanism; that is, cyclic 
creep, which results from a cyclic mean stress and thermal ratcheting, could accel- 
erate failure. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 16, 1976, 

506-21. 
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APPENDIX - COMPARISON OF EQUATION-FACTOR TEMPERATURES 

WITH SINDA TEMPERATURES 

In order to establish that equations (2) to (4) were generally acceptable for predicting 
the radial temperature gradient, temperatures calculated with these equations were com- 
pared with those obtained from SINDA. The thermal maps generated with SINDA for the 
theoretical firing cycle conditions listed in table I were used for comparison. As men- 
tioned previously, equations (2) to (4) provide radial temperature distributions only. 
Therefore, the temperatures calculated with these equations for each row of structural 
elements were considered as average temperatures for the appropriate rows. Element 
temperatures generated with SINDA were averaged for each row and compared with those 
obtained with the equations. 

Figures 12 to 14 show comparisons of element row average temperatures obtained 
with SINDA and those calculated with the equations. The case numbers on each curve 
correspond to a cycle condition listed in table I. Notice that data for two liner materials 
are contained in this table. The thermal properties of these materials are the same so 
that the thermal maps generated for each case are applicable to thrust chambers with 
either liner material. The curves in figures 12 and 13 represent the locus of values ob- 
tained from equation (2), and the curves in figure 14 represent the locus of values ob- 
tained from equations (3) and (4). The symbols associated with each curve are average 
temperatures that were calculated for each row of the appropriate SINDA thermal map 
for each cycle condition. It can be seen that the curves represent the average radial 
temperature distribution quite well. Figure 15 is a comparison of SINDA generated el- 
ement temperatures and those calculated with the equations and their appropriate weight- 
ing factors for the cycle conditions of case 5. The solid line connects the SINDA temper- 
ature for each element, and the dashed line connects those calculated by the equation- 
factor method. The equation -factor curve for the start transient was calculated with 
equation (2) and the weighting factors of table II . The same curve for the shutdown tran- 
sient was calculated with equations (3) and (4) and the weighting factors of table m. Both 
sets of curves compare reasonably well. A structural analysis using the start transient 
equation factor temperatures showed that the strain error of the maximum strained el- 
ement was approximately 1 percent when compared with the maximum strained element 
obtained with SINDA temperatures. Since the shutdown transient curves show better cor- 
relation than the start transient curves, the strain error for these thermal conditions 
should be less. It was concluded then that the equation factor method was satisfactory 
for generating the thermal maps required for the strain and life curves based on the hot- 
gas wall temperature and the temperature difference AT. 
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TABLE I. - THEORETICAL FIRING CYCLE CONDITIONS FOR 
SINDA THERMAL CASES USED TO EVALUATE THE 
APPLICABILITY OF TEMPERATURE EQUATIONS 


[Chamber pressure 4.14 MN/m 2 (600 lb/in. 2 ). ] 


Case 

Liner 

material 

Start 

type 

Cycle 

period, 

sec 

Coolant flowrate 

kg/sec 

lbm/sec 

1 

OFHC a copper 

Slow ramp 

8.0 

0.94 

2.08 

2 

OFHC copper 

Fast ramp 

8.0 

.94 

2.08 

3 

Zirconium copper 

Slow ramp 

8.0 

.61 

1.34 

4 

Zirconium copper 

Fast ramp 

7. 5 

.61 

1.34 

5 

Zirconium copper 

Fast ramp 

1.6 

.61 

1.34 


Oxygen free, high conductivity. 



TABLE H. - MULTIPLYING FACTORS FOR ELEMENT TEM- 


PERATURES AT PEAK AT OF START TRANSIENT 


T i= KT av; 


^_ K 1 +K 2 +K 3+ K 4+ K 5 


Element 

Case 1 
K 1 

csa 

g <N 

3 « 
o 

Case 3 
K 3 

Case 4 
K 4 

Case 5 
K 5 

K 

1 

0.989 

0.991 

1.000 

0.991 

1.000 

0.994 

2 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

3 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

4 

.861 

.869 

.929 

.892 

.918 

.894 

5 

1.039 

1.032 

1.019 

1.029 

1.012 

1.026 

6 

1.103 

1.096 

1.050 

1.083 

1.049 

1.076 

7 

.990 

.990 

.991 

.992 

.992 

.991 

8 

1.009 

1.012 

1.007 

1.008 

1.008 

1.009 

9 

.994 

.994 

.993 

.992 

.994 

.993 

10 

1.006 

1.006 

1.006 

1.004 

1.007 

1.006 

11 

.995 

.992 

.995 

.996 

.995 

.995 

12 

1.005 

1.006 

1.004 

1.004 

1.003 

1.004 

13 

.999 

.999 

.998 

.997 

.997 

.998 

14 

1.001 

1.001 

1.001 

1.003 

1.002 

1.002 

15 

1.046 

1.063 

1.043 

1.053 

1.054 

1.052 

16 

1.034 

1.036 

1.031 

1.038 

1.040 

1.036 

17 

1.004 

1.007 

1.004 

1.004 

1.004 

1.005 

18 

.958 

.956 

.967 

.961 

.960 

.960 

19 

.950 

.937 

.953 

.943 

.941 

.945 

20 

1.039 

1.059 

1.035 

1.042 

1.043 

1.043 

21 

1.033 

1.051 

1.029 

1.035 

1.037 

1.037 

22 

1.002 

1.025 

1.010 

1.012 

1.013 

1.012 

23 

.967 

.964 

.966 

.960 

.959 

.963 

24 

.959 

.940 

.960 

.951 

.950 

.952 

25 

1.028 

1.036 

1.027 

1.032 

1.032 

1.031 

26 

1.023 

1.029 

1.022 

1.026 

1.027 

1.025 

27 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

28 

.978 

.971 

.979 

.975 

.974 

.975 

29 

.972 

.964 

.974 

.968 

.967 

.969 

30 

1.025 

1.032 

1.024 

1.028 

1.024 

1.028 

31 

1.011 

1.014 

1.011 

1.013 

1.031 

1.012 

32 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

33 

.986 

.982 

.987 

.984 

.984 

.985 

34 

.977 

.971 

.979 

.975 

.973 

.983 




TABLE m. - ELEMENT TEMPERATURE CONSTANTS 


AT PEAK AT OF SHUTDOWN TRANSIENT 


[ Ti =CT av or T ay + Cl .] 


Element 

C 

a C x - °R 

1 

1.000 


2 

1.000 

-- 

3 

1.000 

-- 

4 

1.013 

— 

< 5 

.993 

-- 

6 

.991 

-- 

7 


-1 

8 


1 

9 


-1 

10 


1 

11 


-1 

12 


1 

13 


-1 

14 


1 

15 - 34 

1.000 

— 


a Cj represents the variation in temperature of el- 
ements 7 to 14 from the average temperature cal- 
culated with eq. (4). 


TABLE IV. - MAXIMUM STRAIN RANGE VALUES 


Firing cycle case 
from table I 

Total strain range 
with residual 
strains, 

% 

Strain range a 

% 

Start 

Shutdown 

Total 

1 

2.22 

2.08 

0.23 

2.31 

2 

1.95 

1.70 

.20 

1.90 

3 

2.35 

2.22 

.23 

2.45 

4 

2.44 

2.20 

.29 

2.49 

5 

2.31 

2.21 

.21 

2.42 


a From figs. 8 and 9. 







Figure 1. - Fatigue life test thrust chamber. 
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4.03 

1.585 

63 

El 

4.03 
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114 

3.78 
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165 
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1.490 

197 

354 

3.78 

1.490 

209 

376 
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1.430 

345 

620 

3.63 

1.430 
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716 
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Strain range, percent 


^Temperature difference, AT 




Figure 5. - Schematic representation of relation of three cycle 
parameters during two different firing cycles. 



= total strain during start 
= plastic strain during start 
€^ ot "total strain during shutdown 
Figure 6. - Theoretical stress-strain hysteresis loop. 



Figure 7. - Typical low-cycle fatigue life of half hard Zr-Cu. 
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Maximum total strain, percent 



0 -100 -200 -300 -400 -500 

Temperature difference, AT, °F 

Figure 9. - Effect of hot-gas wall temperature and AT on total strain 
of maximum strained element at time of peak AT during shutdown 
transient. 
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Cycles to fracture 













Temperature, 



Temperature 





1.30 1.35 1.40 1.45 1.50 1.55 1.60 1.65 


Radius, in. 

Figure 14. - Comparison of element row average temperatures obtained with SINDA and those 
calculated with equations (2) and (3! for the shutdown transients of cases 1 to 5. 
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Temperature, 



Structural element numbers 

Figure 15. - Comparison of individual element temperatures obtained with SINDA and those calculated by equation 
factor method for case 5. 
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